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LAGEOS is an accurately-tracked, dense spherical satellite covered with 426 retro-reflectors. The 
tracking accuracy is such as to yield a medium term (years to decades) inertial reference frame 
determined via relatively inexpensive observations. This frame is used as an adjunct to the more 
difficult and data intensive VLBI absolute frame measurements. There is a substantial secular 
precession of the satellite's line of nodes consistent with the classical, Newtonian precession due to 
the non-sphericity of the earth. Ciufolini has suggested the launch of an identical satellite (LAGEOS- 
3) into an orbit supplementary to that of LAGEOS-1: LAGEOS-3 would then experience an equal 
and opposite classical precession to that of LAGEOS-1. Besides providing a more accurate real-time 
measurement of the earth's length of day and polar wobble, this paired-satellite experiment would 
provide the first direct measurement of the general relativistic frame-dragging effect. Of the five 
dominant error sources in this experiment, the largest one involves surface forces on the satellite, 
and their consequent impact on the orbital nodal precession. The surface forces are a function of 
the spin dynamics of the satellite. Consequently, we undertake here a theoretical effort to model 
the spin ndynamics of LAGEOS. In this paper we present our preliminary results. 



I. THE LAGEOS-3 MISSION. 



The Laser Geodynamic Satellite Experiment (LAGEOS-3) is a joint USAF, NASA, and ASI proposed program 
to measure, for the first time, a quasi-stationary property of the earth - its gravitational magnetic dipole moment 
(gravitomagnetism) as predicted by Einstein's theory of general relativity. This gravitomagnetic field causes local 
inertial frames to be dragged around with the earth at a rate proportional to the angular momentum of the earth, 
and inversely proportional to the cube of the distance from the center of the earth. Thus the line of nodes of the 
orbital plane of LAGEOS-3 precesses eastward at 32 mas/yr. Although in this example the frame dragging effect is 
small compared to the torque on the orbital plane due to the oblateness of the earth, it is an essential ingredient in 
the dynamics of accretion disks, binary systems, and other astrophysical phenomena JlJ]. 

Today, almost eighty years after Einstein introduced his geometric theory of gravity, we have just begun to measure 
to verify - his gravitation theory. Of no less stature than the "tide producing" —M/r 2 "electric component" of 
gravity is the inertial- frame defining "magnetic component" of gravitation —J/r 3 . To see this force in action: first, 
inject a satellite into a polar orbit about an earth-like mass idealized as not spinning with respect to the distant 
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quasars. The satellite will remain in orbit in a continuous acceleration towards the center-of-mass of the attracting 
body under the influence of the Newtonian 1/r 2 force, and its orbital plane will remain fixed in orientation with 
respect to distant quasars. Second, spin the central body, giving it angular momentum, and follow the trajectory of 
the satellite. Its orbital plane will experience a torque along the central body's rotation axis. The orbital plane will 
undergo a precessional motion in the direction of the central body's rotation. The mass in motion of the central body, 
or "mass current" , produces a dipole gravitational field - the gravitomagnetic field. In the case of a satellite orbiting 
at two earth radii, the orbital plane will precess about the body axis of the earth at approximately 32 mas/yr. This 
is the Lense-Thirring effect Q . 

The Lense-Thirring force has never been directly measured. A measurement of this gravitomagnetic force can 
be compared to the pioneering work of Michael Faraday on the measurement of the magnetic force between two 
current-carrying wires. However, the laboratory setting for this gravity measurement will be the 4-dimensional curved 
spacetime (approximately Kerr) geometry enveloping the earth. The idea behind the LAGEOS gravity measurement 
is simple. Whereas the Everitt-Fairbanks experiment (Gravity Probe-B) proposes putting a gyroscope into polar 
orbit ||], t ne Ciufolini LAGEOS-3/LAGEOS-l experiment Q proposes the use of the orbital planes themselves as a 
gyroscope. 

In 1976 NASA launched the LAGEOS-1 satellite, a totally passive 60 cm diameter ball of aluminum with 426 
retro-reflecting mirrors embedded in its surface. (There are numerous globally-located laser tracking stations to 
observe LAGEOS-type satellites.) LAGEOS-1 was injected into a two earth-radii circular orbit at an inclination of 
110 deg. Due to the oblateness of the earth, the orbital plane rotates at a rate of 126 deg /yr. This torquing can only 
be modeled to 450 mas/yr - which is not accurate enough to measure the 32 mas/yr gravitomagnetic force. The 
idea of Ciufolini is to launch another LAGEOS satellite (LAGEOS-3) into an orbit identical to that of LAGEOS-1, 
except that its inclination is supplementary (70 deg = 180 deg— 110 deg). This proposed orbital plane will rotate 
in the opposite direction, i.e., —126 deg /yr. The intersection of the two (LAGEOS-1, LAGEOS-3) orbital planes 
will sweep out a "tandem-generated gyro plane" (Fig. [J). The utilization of two satellites cancels out many of the 
large precessional effects due to mass eccentricities of the earth, providing a plane inertial enough for a measurement, 
accurate to five percent (or better) , of the "magnetic component" of gravity. 




FIG. 1. The tandem-generated gyro plane and the LAGEOS-3 satellite experiment. Both the Gravity Probe B 
(GP-B) experiment at Stanford and the LAGEOS-3 experiment share the feature of generating an effective polar orbit so as to 
cancel the torque on the orbital plane arising from the oblateness of the earth. In the case of the GP-B experiment, a drag- free 
superconducting gyroscope is launched into a polar orbit. In the case of the LAGEOS-3 experiment, a satellite is launched into 
an inclination that is supplementary (70 = 180 — 110 deg) to the LAGEOS-1 satellite. As the orbital planes precess in equal 
and opposite directions, the line of intersection of these two supplementary planes evolves so as to sweep out a polar plane. 
Einstein's theory of general relativity predicts that this tandem-generated gyroplane will precess at a rate of 32 mas/yr. In 
the GP-B case, the gyroscope is ~ cm in diameter, while it is 12, 271.79 km in diameter for LAGEOS. 

II. WHY SPIN DYNAMICS? 

The success of the LAGEOS experiment hinges upon the detection of a 32 mas/yr eastward drift of the line of 
nodes of the two satellites. A strong effort is now underway to model the orbital and spin dynamics of these satellites, 
and to make an assessment of the uncertainties these will add to the desired measurements. Encapsulated in Table ffl 
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are the five major classes of errors in this experiment: (1) geopotential (other than even zonals) and tides, (2) earth 
radiation pressure, (3) uncertainty in other relativistic effects, (4) earth and solar-induced thermal forces, and (5) 
even zonal geopotentials (per 0.1 deg inclination injection error). Due to the recent GEM-T1 improvements in the 
earth's zonal harmonics, the errors due to solid-earth tides have been significantly reduced, and are now potentially 
smaller than those due to surface effects fa] . 

These surface effects (e.g., Yarkowsky thermal drag, neutral and charged particle drag) cause a change in the nodal 
precession of the satellite, thus contributing to what will be potentially the largest source of error in the LAGEOS-3 
measurement of the gravitomagnetic effect. In modeling these surface forces, in particular the Yarkowsky thermal 
drag (also referred to as the Rubincam effect) due to the satellite's differential heating and delayed reradiation, the 
behavior of the spin vector of the satellite is of crucial importance. That the uncertainties induced by the surface 
forces on LAGEOS arc on the order of 4%, out of a 6% experiment, makes such a theoretical modeling of paramount 
importance. With this in mind, we undertake in this paper a theoretical model of the spin dynamics of LAGEOS-type 
satellites, and compare our results with observational data. 

We must emphasize that we are looking for qualitative, not quantitative, results. We do not intend to predict the 
exact magnitude and direction of the spin vector at any particular instant. Instead, we desire a rough estimate of 
its magnitude, and whether its behavior is predictable or chaotic. We seek to answer questions of the sort: "In the 
asymptotic limit, does the satellite tidally lock, tumble, or have some other behavior?" 

Previous studies of the spin dynamics of the LAGEOS-1 satellite [Q were valid only for spin rates much larger than 
the orbital frequency. Today, however, the spin of the satellite, decaying with a three-year time constant, is rapidly 
approaching the orbital period. We therefore require an analysis in the low-frequency regime. Unlike the earlier 
analysis which dealt with orbit-averaged quantities, we will solve numerically the full set of dynamical equations. As 
we will show, while the low frequency regime exhibits complex behavior, the asymptotic state of this forced, damped 
system appears to be tidally locked. In this paper we will examine the spin-orbit resonance phase and discuss the 
asymptotic state of the spin of the satellite. It is the goal of this research to provide an optimum strategy for the 
measurement of the spin dynamics of the LAGEOS satellites in support of the proposed gravity measurement. In 
addition, we use our theoretical model to propose an optimum orbital injection procedure for LAGEOS-3: it is our 
opinion that the LAGEOS-3 satellite should be injected into the orbital plane with as large a spin rate as possible. 
Our results provide the first analysis of the asymptotic spin dynamics of the LAGEOS satellite. Previous calculations 
were unable to analyze the spin-orbit resonances of the satellite, nor its asymptotic behavior, which will play a crucial 
role in the experiment. We demonstrate in this paper that the LAGEOS-1 satellite will be sufficiently predictable to 
support the gravitomagnetic measurement. 



Geopotential and tides 


2% 


Earth radiation pressure 


1% 


Uncertainty in other relativistic effects 


1% 


Thermal forces 


1-3% 


Even zonal geopotential 


<1% 


Random and stochastic errors 


2% 


RSS error 


3-4% 



TABLE I. The latest error budget for the LAGEOS-3 Lense-Thirring experiment. Errors in geopotential and 
tides reflect improvement over GEM-T1. Uncertainty in the thermal forces depends mainly on knowledge of spin-axis motion. 
The even zonal geopotential error calculation assumes less than 0.03 deg inclination injection error. The random and stochastic 
error estimate can accommodate seasonal variations in low degree spherical harmonics of the geopotential. 
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III. MODELING THE SPIN 



It is rather interesting that after 36 years the spin dynamics of passive satellites is once again important to the 
field of astrodynamics. In 1957 Vinti S analyzed the spin dynamics of a non-ferromagnetic spherical satellite in the 
earth's magnetic field, which was then of critical importance to the alignment of antennas. Today we perform the 
same analysis on a slightly oblate satellite, of critical importance to the first measurement of the magnetic component 
of gravity, as predicted by Einstein. 

There are many factors to consider when analyzing the spin dynamics of an oblate, metallic satellite orbiting in the 
gravitational and magnetic fields of the earth. The most prominent effect is the torqueing due to the gravitational field 
of the earth. This arises from the oblateness of the satellite, with the oblateness of the earth producing a negligible 
contribution that can be added to our calculations as an adiabatic correction. If the satellite's (bulging) equatorial 
plane lies in the plane of orbit, no such torques are possible. However, when the satellite is not placed exactly in such 
a position, gravitational torques will arise. In an effort to model these torques, we consider the situation of an oblate 
spheroid in orbit around a point mass. As was done by Bertotti and less, we parallel the development in Goldstein 



Bertotti and Iess's analysis of the effects of gravity on the spin dynamics of the oblate satellite, which resulted in 
predicting a chaotic spin dynamics of arbitrarily large amplitude in the obliquity at late times, is not appropriate for 
small rates of spin. Their prediction is based on the "Hipparcos" formula for the rate of precession uj p of an oblate 
satellite in the inhomogeneous gravity field of the earth 

o; p = ^A^°cos0, (1) 

where ujq is the orbital angular velocity, ui^ is the satellite spin rate, is the obliquity angle of the satellite (the angle 
between Q and the normal to the orbital plane) and 

A = ^ (2) 

is the satellite's oblateness. Here I3 and I\ — I2 are the principal moments of inertia (the principal direction 
corresponding to ^3 is assumed to be that of Q and coincides, by assumption, with the symmetry axis of the oblate 
satellite). 

It is argued by Bertotti and less that, since u p oc l/w 3 , the gravitational precession in the asymptotic limit of small 
lo becomes very fast and may make the spin dynamics chaotic. This conclusion is based upon a misunderstanding. We 
have shown, via a careful analysis of assumptions underlying the "Hipparcos" formula, that even in the approximation 
commonly used in deriving the formula (averaging of the gravitational potential over the satellite orbit, dipole cutoff 
of the multipole decomposition, etc.) Eqn. ([!]) can only be used when 

6^gA cos 2 9 w 2 

, ,2 ~ , ,2 ^ x ' \ 6 ) 
uj 3 uj 3 

i.e., when the spin rate of the satellite is much greater than its orbital angular velocity. The latter restriction is easy 
to overcome, and the corrected equation for u> p is 



This equation imposes a bound on 



/3 

\uj p \ < uj J -A, (5) 

which makes it clear that the rate of precession cannot grow to cause chaoticity of the satellite spin dynamics. A 
subsequent qualitative investigation jlCj] by Chris Fuchs has shown that when magnetic forces are included in the 
picture the precession rate remains bounded, and should be much smaller than luq. Another conclusion reached in 
the course of our analysis has been that, when both gravitational and magnetic forces are taken into account, the 
nutation, although bounded in its amplitude, does not disappear completely even in the asymptotic limit. The value 
of these results, however, is limited by the fact that they do not yield the exact bounds, nor do they provide any 
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information on the time scale to reach the asymptotic limit. However, they do clearly show that chaoticity of the 
LAGEOS spin dynamics caused by an unbounded growth of gravitationally-induced precession cannot occur. 

Another important factor governing the evolution of the spin vector is the interaction of the metallic core of the 
satellite with the magnetic field of the earth. The LAGEOS satellite Jll| is composed of two aluminum hemispheres 
bolted together, with a brass cylindrical core along its body axis (the original axis of spin) (Fig. §. The spinning 
of this metallic object in the magnetic dipole field of the earth (and the motion through that field) will cause eddy 
currents within the satellite, which will in turn cause dissipation through Joule heating and a slowdown of the spin, 
and furthermore will cause torques on the spin vector. These torques can be understood as the interaction of the 
magnetic dipole, caused by the induced eddy currents, with the magnetic field of the earth. In modeling this effect, 
we have treated a uniform, spherical object in the orbit of a perfect magnetic dipole. Here wc are concerned with 
a simple qualitative analysis of the spin dynamics of LAGEOS, as the complexity of the true satellite geometry 
prohibits us from a precise model of the eddy current distribution. Our theoretical model permits us to analyze 
arbitrary inclinations of the orbital plane to the earth's magnetic dipole axis. The results presented in this paper 
correspond to a retrograde I = 109.86 deg orbit. 

IS-1/94-853 4/93 




FIG. 2. A schematic representation of the LAGEOS assembly. The two hemispheres form a 60 cm diameter sphere 
and are made of 6061 aluminum. There are 426 cube corner reflectors recessed into the satellite's surface. To obtain the desired 
surface-to-mass ratio, the inside of the aluminum hemispheres were hollowed out to make room for a brass cylinder of 31.76 cm 
diameter and 26.70 cm height. The LAGEOS-1 satellite was launched June 7, 1976 into an inclination of 109.8585 deg and 
semimajor orbit of 12271.790 km. (This figure is adapted from Ref. [juj.) 

The problem of a spinning metallic sphere in a constant magnetic field has been treated by Landau and Lifshitz [fl2"| , 
and we avail ourselves of their results. For our purposes, we ignore the torques caused by the changing magnetic field 
due to the orbit (as opposed to spin) of the satellite. These torques can be shown to be negligible until asymptotically 
late stages of motion, and have no qualitative effects upon the dynamics. 

The source of many of the difficulties in doing analyses of such orbiting, spinning bodies lies in the involved 
coordinate systems needed to describe their motion. Thus, it is important at this point to give a brief description of 
the coordinates we will use in this paper. In our analysis of the spin dynamics of the LAGEOS satellite we found it 
convenient to introduce the following four coordinate systems: 

1. The orbit- centered inertial frame ( OCI) {xi, y±, z\}. Here the zi-axis is oriented along the normal to the orbital 
plane of the satellite (~ 70.142 deg coinclination) . The Xi-axis is defined to be the intersection of the orbital plane 
and the earth's equatorial plane, and the origin is the center of mass of the earth. We have assumed here that this 
frame is inertial, and have not included the secular drag of the line-of-nodes of the orbital plane due to the oblateness 
of the earth. This ~ 126 deg jyr precession can be included at the end of our analysis as an adiabatic correction, and 
again has no qualitative effects upon the dynamics. 

2. The earth- centered inertial (ECI) reference frame \xi-,yi-, z-{\ . Here the ^2-axis is aligned with the body axis of 
the earth. The X2-axis lies in the earth's equatorial plane at zero degrees longitude, and the origin coincides with the 
center of the earth. 

3. The body frame (non-inertial) {xb, yb, Zb}- The origin is at the center of the satellite, and the axes correspond to 
a set of principal axes. The satellite is assumed to be a slightly oblate (~ 3.27%) spheroid of brass, and the z&-axis is 
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aligned along its body axis. The Xb and yb axes are an arbitrary fixed set of orthogonal axes spanning the equatorial 
plane of the satellite. In our calculations, the body axis is related to the orbit-centered frame through the three Euler 
angles 9, (p, and ip. The nutation angle 9 is the angle between Zb and z±, while the angle of precession <f> is the angle 
between the xi-axis and the line of nodes (the intersection of the orbital plane and the equatorial plane of the satellite 
{xb-yb plane)). The spin angle ip is the angle between the line of nodes and the a^-axis. 

4. The Landau- Lifshitz (non-inertial) coordinate system {x ll ,y 11 , z 11 }. The z H -axis is aligned along the instantaneous 
angular momentum vector of the satellite (tD). The x"-axis is picked so that the instantaneous magnetic field (B) at 
the satellite lies in the x 11 — z 11 plane. Note that z 11 need not be aligned with the body axis of the satellite, and in 
fact, during the asymptotic behavior of the satellite they are vastly different. In particular, the angle (£) between the 
symmetry axis of the satellite (zb) and the instantaneous angular momentum vector (zu) can be expressed in terms 
of the three Euler angles as 

sine=fl + |^±itV 1/2 (6) 
^ II (0 2 sin 2 6 + 02 ) J 1 ; 

In the early stages of the LAGEOS mission when ip » 9 and ip ^ <t> this angle is rather small (£ w (Ii/I^fae/u^) w 
6.3516 x 1CT 10 ). 

A proper interpretation of the results of the numerical simulation requires a careful distinction between to p and u)^ 
(where <fi is the Euler angle of the body frame). The two coincide only under the assumption that the oj^ component 
of the satellite's angular velocity is responsible for all of the satellite's energy, or, to put it differently, the angular 
momentum of the satellite is directed along the body axis of the satellite. As the satellite's spinning motion slows 
down, this last assumption is violated. uj p represents only a part of ojq, the other part being caused by the "tidal 
locking" effect. As we shall see, this is exactly what happens in the asymptotic phase of the satellite's motion. 
Nevertheless, both uo p and lo^ remain bounded. 



IV. SPIN DYNAMICS OF LAGEOS: THE EQUATIONS. 

The spin dynamics are determined by Euler's equations 

/iWi - u 2 Us{h ~ h) = Ni, 

h^2 ~ U3Ui(h -h) = N 2 , (7) 
I 3 uj 3 = N 3 , 

where w\, u>2, ^3 are the components of the satellite's angular velocity in the body frame, I\ = I2, I3 are the principal 
moments of inertia, and Ni, N 2 , N 3 are the components of the torques along the satellite's principal axes. After 
substituting the expressions for w\, 0J2, ^3 m terms of the Euler angles 

oj\= 4> sin 9 sin ip + cosip, 

uj 2 = <j>sin9cosip — 9 simp, (8) 
= <j)COS0 + ip , 



the Euler equations become 



•a a , Nicosip- N 2 smtp 

V — Vfree H ; , (.» J 

J-l 

9 f ree = (— —J (fP cos 9 sin 9 — —ipcpsm9, (10) 

V h J h 

2 2 ^ sin ip + N 2 cos ip 

<p = (pfree-\ T . Q , (H) 

.. ( h-2h \ cose-. i 3 eip 

<Pfree=[ ; ^^^^+7"^^' ( 12 ) 

\ h J sin I\ sin 9 

N 3 Ni sin ip + N 2 cos ip cos 9 . . 

* = *tr~ + T Ti — 9 , (13) 

tpfree = - 7 — — TT^Cp + ~ ~T - ~T ■ (14) 

\ Ji / sinO smt) h sin 6 
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The torque components N±, N2, A3 are due to gravitational and magnetic forces acting on the satellite 

N i = N$ s) +N$ m \ i = l,2,3. (15) 
Gravitational torques in the body frame are given by 

i„) dV sini/> dV cos 9 simp dV 

N\ 9) =- C0S tP— - -^jrr + 

ou sin 9 dtp sm0 dtp 



( ) dV cosip dV cos 9 cos ib dV 

06 smfc> d(f> sinf dip 



N- 



(9) 



dV 



3 dip' 

Using the standard dipole approximation, the gravitational potential V is 

V = -1), (IT, 

where 7 is the direction cosine between (1) the radial vector from the satellite center of mass to the center of the 
earth, and (2) the symmetry axis of the satellite. It is related to Euler's angles via 

7 = sin 9 sin (77 — 770 — <p), (18) 

where 77 gives the angular position of the center of mass of the satellite in its orbit about the earth, and 770 is an 
arbitrary starting position. The potential is given by 

V = GM 2 3 R 3 h) (3 «in 2 sin 2 ^ - m - <f>) - l) , (19) 



and thus 



8V 3GM(I 3 -h) 2 

— — sin 9 cos 9 sm (77 — 770 — 



89 R 3 
1 dV _ 3GM(I 3 - h 



■ shi9sm(rj — 770 — <p) cos(?i — 770 — (p), (20) 



sinfl dcp i? 3 

%-* 

Equations (10) and (14) lead to the following final expressions for the components of the gravitational torque in the 
body frame 

(s ) 3GAf(/ 3 -7i) . . f -cos 9 cosip sin (77- 770 

N i = ^ sin0 sm 77 - 770 - <p) < . . , , 

f g ) 3GM(h~h) . . f cos 6» sin ^ sin(?7 - 770 - . 

^2 = ^5 sm sm(77 - 770 - <p) i . , / , 

J J?3 u r/ [ + cost/; cos(t7 — 770 — (p 

A 3 (s) = 0. (21) 

As for the magnetic effects, we are interested in the torque components acting on a conducting ball of radius a 
spinning with angular velocity Co in an external magnetic field B. We take our expressions from Landau and Lifshitz 
[ p2[ , noting that their results are for what we have dubbed the Landau-Lifshitz frame (as described above), and for 
our purposes need to be transformed to the body frame. 

As already described, the Landau-Lifshitz frame is determined by the vectors Co and B. If we introduce an arbitrary 
set of rectangular coordinates x, y, z, in this frame Co and B are represented as 

Cd = LO x X + LO y y + LO z Z = (tO x ,LOy,LO z ), (22) 

B = B x x + B y y + B z z = (B x ,B y ,B z ), (23) 

where hats (x) denote vectors normalized to unity. The transition between this arbitrary frame and the Landau- 
Lifshitz frame is given by 
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o)|wx B\ 



V 



[(SxB)xlu]. 
u)\ujxB\ 



uj\SxB\ 
("XB)» 



|c3xB| 



/ 




(24) 



where lu = \uj\. 

The components of the magnetic torque in the Landau-Lifshitz frame are given by 



Va"B%B l z l , 
~Va'B"B u 



N u 



N l J = -Va"(B») 



Z 1 



(25) 



where V = 47ra 3 /3 is the volume of the ball, By, B l z l are components of the magnetic field in the Landau-Lifshitz 
frame, and the real and complex parts of the coefficient of magnetization are 



a = 



_3_ 

"8tt 



a" = - 



98 2 



IGira 2 



36 sinh (2f) - sin (2f ) " 
2a cosh (2f) - cos (2f) 

i o sinh (2f ) + sin (2f ) 
6cosh(2§) -cos(2f) 



with 



\J2-kouj 

Here c is the speed of light and a is the specific conductivity of the material forming the ball. At small values of lo 
(5 » a), a' and a" can be approximated by the expressions 



105 c 4 
a 2 auj 



a 



10c 2 



(26) 
(27) 



Components of the magnetic field can be evaluated in the ECI frame using the dipole approximation. In the 
spherical coordinate representation of the ECI frame, the magnetic field components are 



B r 



2M 



COS 01, 



B 6l = --^ sin 0i, 



(28) 



Ba 



o, 



where M is the magnetic dipole moment of the earth. In the rectangular coordinate representation of the ECI frame 
(the index 1 is used everywhere for quantities in ECI) 



B x 

By 

B 7 



B ri sin 0i cos 4>\ + Bg 1 cos 0i cos </>i , 
B ri sin 0i sin 0i + Bg 1 cos0i sin^i, 
B ri cos 0i — Bg 1 sin 0i. 



(29) 



Transforming these to the OCI frame: 



B x 
B„ 



{B r 
{B r 



Bg 1 COS 01 ) COS (f)\, 

Bg 1 cos 0i ) sin 0i cos£ + (B ri cos0i 



9 1 sin0i), sin£ 

B Z2 = — (B ri sin 0i + Bg 1 cos 0i ) sin <j>i sin £ + (B ri cos 0i — Bg 1 sin 0i ) cos £, 



(30) 



where £ is the colatitude angle of the orbit with respect to magnetic north, and the index 2 is used for quantities in 
the OCI. The satellite's semimajor axis is R — 1.227179 x 10 9 cm, and for the earth's magnetic dipole moment we 
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used M = 7.9 x 10 25 G — cm 3 . Although we do not average the magnetic field in our simulations, it was useful for 
back of the envelope calculations to note that the the magnetic field, averaged over one orbit, is \B\ = 0.155814 G, 
with components, 

B X2 = G, 

B y2 = -0.0388558 G, (31) 
B Z2 = -0.150891 G. 

The angles 0\ and 0i are dependent upon the satellite's position in its orbit. The satellite's coordinates in the OCI 
frame are 

x 2 = rcos(?7 - t/o), 

y 2 = rsin(ry - no), (32) 

Z2 =0, 



where 



In ECI we have 



T, 



2vr , s 

t. (33) 



orbit 



xi = Rcos(r/ - no), 

yi = R sin(ri - no) cos f , (34) 
Z\ = — i?sin(7y — no) sin£. 



Hence, 



sin 0i = \J cos 2 (n — no) + sin 2 (?7 — 770) cos 2 £, (35) 

cos 9\ = — sin(?7 — 770) sin£, (36) 

COs(r; - no) ,o 7 ^ 

cos 0i = — , (37) 

sinfi 

. sin(?7-7/o)cos$ 

sin 0i = — . 38 

The three Euler equations (Eq. |l4]), with the magnetic and gravitational torques included (properly transformed 
to the body frame), give us the vehicle to analyze qualitatively the spin dynamics of the LAGEOS-1 satellite. We 
present our results from the numerical integration of these equations in the following two sections. 



V. INITIAL- VALUE DATA. 

We solved the Euler equations (Eq. |l4|) using a fourth-order Bulirsch-Stoer algorithm with adaptive time-stepping 
jn| . The equations were integrated for 3 x 10 9 seconds, as we wanted to (1) reproduce the experimentally-measured 
spin rates (17 yrs following launch), (2) examine the spin-orbit resonance (~ 27 yrs after launch), and (3) reveal the 
asymptotics of the spin dynamics (~ 79 -I- yrs after launch). The experimentally- measured exponential decrease in 
the spin rate imposed a constraint on our theoretical model, linking the "effective" radius of the satellite (a) with the 
satellite's "effective" conductivity (a): 



a a 



5 - 1.195 52 x 10 24 cm 5 /s. (39) 



The satellite was modeled as a 25.55 cm radius spheroid of brass (a = 1.098 x 10 17 s _1 ). LAGEOS Ps moment of 
inertia about the body axis is I Zb = 1.314 x 10 s g — cm 2 , while the moment of inertia perpendicular to the body axis 
is I Xb = I yb = 1.271 x 10 s g — cm 2 , corresponding to an approximately 3.38% deviation from sphericity. 

From the experimental data, we observed a deviation from pure exponential damping of the spin of the satellite 
at early times. This is presumably due to the satellite's transition from magnetic opaqueness to transparency in the 







course of its spin damping. A rapidly rotating conductor, with angular velocity oj and conductivity a, will have an 
associated magnetic skin depth 



5 = c/^/2-Kau). 



(40) 



This skin depth starts out considerably smaller than the satellite, but as the satellite's spin is damped down, the skin 
depth becomes much larger than the satellite's diameter. Although this transition effect will be more pronounced 
for our idealized spherical brass model satellite than for LAGEOS (with its additional surface structure), it helps 
account for the structure in the experimental spin-down data (for which it would have been convenient to have access 
to the associated error bars). It is for this reason that we began our calculations within the transparent phase, at 
t = 92772864 s. We integrated forward in time to t = 3 x 10 9 s, and backwards to t = s. At t 0} the skin depth 
of our satellite model, given the experimentally measured spin rate of ui n = 4.36332 and a conductivity for brass 
(a = 1.098 x 10 17 s _1 ), is <~ 17.3 cm. We started the integration with the satellite in a retrograde circular orbit, at 
a radius of 1.227179 x 10 9 cm and an inclination of 109.859 deg. At t 0} the satellite was located over the equator 
and positioned along the x-i — x\ axis = 0). We inserted the body axis of the satellite into the orbital plane with 
the angular momentum vector parallel to the OCI yi-axis. The initial spin of the satellite was determined by the 
experimental measurements, and was set to 4.36332 rad/sec. The initial conditions were 



In running this initial data backward in time from t = 927728640 s to launch at t — 0, we recovered the following 
angles and angular velocities (measured in s _1 ): 



The magnetic field in our simulation was assumed to be a perfect dipole field, of moment —M = 7.9 x 10 25 G/cm 3 
and aligned along the z 2 -axis of the ECI frame. 



tp(t ) = 0, 0(i o ) = 1.5707963, <j>{t ) = 0, 
6(t Q ) = 0, <p(t ) = 0, ip(t ) = 4.36332 s~ 



-l 



(41) 



6{t = 0) = 1.304029, 0(0) = -1.327836, V(0) = -0.5542969 

9(0) = 6.792913 x 10~ 10 , <j>{0) = 3.806309 x 10~ 9 , ip(0) = 10.77009. 



(42) 
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FIG. 3. The spin dynamics of LAGEOS: Numerical results. We present here results from our numerical simulation 
of the dynamic evolution of the LAGEOS satellite (here modeled as a slightly oblate spheroid of brass orbiting an earth-like 
mass idealized as perfectly spherical and endowed with a polar-oriented dipole magnetic field, top left figure). The evolution 
of the satellite's angular momentum reveals three unique phases, as demonstrated in the plot of the component of angular 
momentum orthogonal to the orbital plane (top right). The first phase is characterized by an exponential decrease in spin of 
the satellite, with negligible nutation and precession. The decay in spin for a 25.55 cm radius sphere of brass agrees with the 
data (bottom left). In the spin-orbit resonance phase the angular velocity decays to a value comparable to the orbital angular 
velocity. The nutation and precession increase, and the angular momentum vector lifts off the orbital plane and settles down 
orthogonal to the orbital plane (top right). The third phase (asymptotic) is characterized by the tidal locking of the satellite. 

VI. RESULTS. 

We have identified three distinct phases in the spin dynamics of the LAGEOS satellite (Figs. 3-5), which we shall 
refer to as (1) the Fast Spin Phase, (2) the Spin-Orbit Resonance Phase, and (3) the Asymptotic Phase. 

a) Fast Spin Phase. This first phase, from t = to t ~ 25 yrs, is characterized by an exponential decrease in the 
spin rate (ij>) of the satellite (Fig. |) from ip(t = 0) = 10.7709 s" 1 to ip(t ~ 25 yrs) w .001 a -1 . The body axis of 
the satellite is aligned with the angular momentum of the satellite in this phase; hence, all other quantities (angular 
momentum, kinetic energy, etc.) also decrease exponentially. The nutation of the satellite (the angle between the 
body axis of the satellite and the normal to the orbital plane) undergoes a steady increase from 9q = 90 deg at t = to 
to 9 ~ 126 deg over a period of ~ 6.4 yrs, indicating a nutation angular velocity of Cug sa 20 deg jyr. The nutation 
then settles down to a "pseudo-stable" state at 9 ~ 115 deg over the next 3.2 years and remains at this value (±5%) 
until the onset of the spin-orbit resonance at ~ 25 yrs (Fig. |6j). Finally, within this twenty five year period the satellite 
precesses in a positive sense by 63 revolutions before unwinding when entering the spin-orbit resonance phase. 

b) Spin-Orbit Resonance Phase. The spin dynamics abruptly change at t ~ 27 yrs. This is precisely when the 
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satellite's spin, decreased by magnetic damping, approaches the orbital angular velocity 



/ 3 9 x 10 20 

ow = y (1 . 227179 xl0 9 )3 s l - 0.0004593782 s~\ (43) 

The conductivity of the satelite was chosen to reproduce experimentally observed exponential decay in the spin rate 
of the satellite, 

i>{t)=i>{t a )e-^, (44) 

with (A ~ 2.96 yr time constant). Therefore, the spin-orbit resonance should occur at t sp ~ 30 yrs, which is in 
agreement with the numerical results. The resonance phase is marked by a movement of the angular momentum 
vector to a position orthogonal to the orbital plane, and is furthermore characterized by the beginning of satellite 
wobble (i.e., the point in time when £ of Eq. (H) becomes nonzero and the body axis becomes misaligned with the 
instantaneous angular momentum vector). From this point forward it is more illustrative to examine the dynamics of 
the instantaneous angular velocity and momentum rather than the Euler angles. In addition to the dramatic changes 
in the spin dynamics, this second phase also gives rise to a reversal in the signs of the precessional velocity (ujq) and 
spin (aty). 

c) The Asymptotic Phase. Following the spin-orbit resonance phase, the spin dynamics gradually settled down to 
an asymptotic regime over the course of ~ 50 yrs. Not surprisingly, the satellite becomes tidally locked (think of 
the moon). In particular, the asymptotic value of the total angular velocity is equal to the orbital angular velocity, 
subject to small fluctuations. We note that in this phase the torques induced from the changing magnetic field due to 
orbital motion of the satellite will become important, and should no longer be ignored. However, their addition will 
not significantly change dynamics, as the energies at this point are quite low. 

From our numerical runs, a rough estimate of the asymptotic behavior of the satelite model (modulo phase and a 
finite offset in (ft) is given by 

6{t) ~ 1.57 + 0.12 cos(4.6 x 10~ 4 i), (45) 
(j){t) ~ -4.59 x 10-H, (46) 
i){t) ~ 3.04 + 0.115 cos(4.6 x 10~H), (47) 
Lj e {t) ~ 5.5 x 10~ 5 cos(4.6 x 10 _4 i), (48) 
L04,(t) ~ -4.59 x 10~ 4 + 3.1 x 10~ 3 cos(1.4 x 10~H), (49) 
aty(t) ~ 5.3 x 10 -5 cos(4.6 x 10~ 4 i). (50) 



The asymptotic values of the other relevant parameters (Fig. ||) 



K.E. ~ 13.42er ff s, (51) 

- (MB) ~ 3.2 x 10- s ergs, (52) 

G.E. 0.45375er 5 s to - 0A5ergs, (53) 

E - 12.96er#s, (54) 

L - 5.84 x 10 4 g cm 2 /s, (55) 

4 ~ 0.115, (56) 

L y ~ -0.15 ±0.05, (57) 

L z - -0.9933, (58) 

lj - 0.0004595s" 1 (59) 
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FIG. 4. The behavior of the angular momentum vector and its components throughout the mission. The 

evolution of the angular momentum reveals the transition to an asymptotic phase wherein the angular momentum vector is 
roughly orthogonal to the body axis of the satellite, and is consistent with a precession rate equal to the orbital angular velocity. 
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FIG. 5. Dynamics of the Euler angles. The evolution of the nutation (angle of obliquity, 9) of the model satellite 
(upper left). In the asymptotic limit the angular velocities of nutation and precession average to zero (upper and lower right, 
respectively). The precession rate (w^), on the other hand, locks into the orbital velocity (lower left). This last plot demonstrates 
clearly the dynamics through the spin-orbit resonance phase. 
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FIG. 6. Energy. A plot of the logarithm of the total energy of the satellite as a function of time reveals an asymptotic 
limit of 12.96 ergs consistent with tidal locking. The linear region early on (0 — 9 x 10 s s) reflects the resistive damping phase 
of the induced eddy currents. Asymptotically, the satellite orients itself so as to minimize the gravitational potential energy 
(lower left) while the magnetic energy (~M ■ B) becomes negligible. 



VII. CONCLUSIONS. 



Of the five largest sources of error identified in the LAGEOS-3 experiment, the earth and solar-induced surface forces 
are potentially the most troublesome (Table The anisotropic heating of the satellite, and subsequent reradiation, 
gives a "thermal rocketing" perturbation (referred to as the Rubincam effect or the Yarkowsky thermal drag) which 
tends to degrade the experiment. To model this effect requires, in part, a detailed knowledge of the behavior of the 
angular momentum of the satellite. Toward this end we have derived, and solved numerically, a simplified set of Euler 
equations that evolve the angular momentum vector for a slightly oblate spheroid of brass orbiting an earth- like mass, 
idealized as being a perfect sphere and having a perfect polar-oriented dipole magnetic field. The Euler equations 
included both the tidal gravitational torques and the eddy-current torques, as well as the resistive damping torques, 
as modeled by complex magnetization coefficients. Using this rather simplified model, we have identified three phases 
of the rotational dynamics - a fast spin phase, a spin-orbit resonance phase, and an asymptotic phase (Fig. |3|). We 
have also identified an error in the previously established model of asymptotic spin dynamics |7j. This error has led 
to confusion and, in attempts to reconcile observed data with theoretical predictions, has led others to hypothesize 
erroneous models for the moments of inertia of LAGEOS-1 p3|. 

Our results have led us to formulate four as yet unresolved questions: (1) Can we obtain the asymptotic solution 
analytically, and in so doing can we understand the wobbling or slippage of the Euler angles with respect to the 
relatively stable total instantaneous angular velocity?; (2) Can we understand why the rms fluctuations in the grav- 
itational potential energy cause the asymptotic value of the angular momentum vector to be offset from the orbital 
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plane by ~ 10 degl; (3) Can we understand why the nutation angle (9) drifts initially at a rate of ~ 20 deg/yr and 
reaches a pseudo-stable value of ~ 115 degl; and (4) Can we understand the fluctuations in the spin rate (ip) over the 
first ~ 20 yrs, which do not appear to have been detected experimentally? We are addressing these questions by (1) 
introducing a more realistic model of the satellite and earth into our calculations Jl5| ], and (2) exploring more of phase 
space by way of Poincare sections. The results presented here provide us with clues that must be pieced together to 
reveal the physics behind the complex motion we observe. 

The current spin dynamics model suggests that we launch LAGEOS-3, with as large a spin (ip) as possible, into 
an obliquity of 9 ~ 115 deg; although, due to the qualitative nature of our results, the precise numbers are far from 
being conclusive. 

We are currently working with the Center for Space Research at the University of Texas at Austin to determine 
the impact this revised model of the spin dynamics of LAGEOS will have on the LAGEOS-3 mission (in particular, 
how will the Rubincam effect alter the line of nodes of the orbital plane?). In addition, we are working closely with 
colleagues at the University of Texas and the University of Maryland to reconcile the experimental measurements of the 
spin dynamics of LAGEOS-1 with our theoretical model fl6| . Furthermore, we will propose an optimal experimental 
measurement schedule in support of the proposed LAGEOS-3 mission. 
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